Space-time Phase Transitions in Driven Kinetically Constrained Lattice Models 
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Kinetically constrained models (KCMs) have been used to study and understand the origin of 
glassy dynamics. Despite having trivial thermodynamic properties, their dynamics slows down dra- 
matically at low temperatures while displaying dynamical heterogeneity as seen in glass forming su- 
percooled liquids. This dynamics has its origin in an ergodic-nonergodic first-order phase transition 
between phases of distinct dynamical "activity" . This is a "space-time" transition as it corresponds 
to a singular change in ensembles of trajectories of the dynamics rather than ensembles of config- 
urations. Here we extend these ideas to driven glassy systems by considering KCMs driven into 
non-equilibrium steady states through non-conservative forces. By classifying trajectories through 
their entropy production we prove that driven KCMs also display an analogous first-order space-time 
transition between dynamical phases of finite and vanishing entropy production. We also discuss 
how trajectories with rare values of entropy production can be realized as typical trajectories of a 
mapped system with modified forces. 

PACS numbers: 64.70.Pf, 05.40.-a 



Introduction. - The transformation of a supercooled liq- 
uid into an amorphous solid glass at low temperatures 
has fascinated scientists for decades pQ, yet a compre- 
hensive theory explaining the microscopic origins of this 
"glass transition" is still missing. Several scenarios are 
advocated in the literature. These include: (i) A tran- 
sition to an "ideal" glass state reached at some critical 
value of a thermodynamic parameter such as the tem- 
perature. This transition may be thermodynamic, as in 
the random first-order transition theory [2J, or kinetic, 
as in mode-coupling theory [5]. (ii) An "avoided" glass 
transition, controlled by an unreachable thermodynamic 
critical point pQ. (hi) The glass transition as a purely 
dynamical phenomenon that is not controlled by struc- 
tural changes and not related to a phase transition in a 
thermodynamic sense. This approach jS] considers dy- 
namic heterogeneity [5J, the observation of large spatial 
and temporal fluctuations in "dynamical activity" , to be 
the crucial feature of glass formers and derives most of 
its predictions and insights from the study of kinetically 
constrained models (KCMs) of glasses [TJ. 

Analytical studies of KCMs |Sj and numerical investi- 
gations of atomistic glass forming fluids j9] have shown 
the existence of a true phase transition in ensembles of 
trajectories. This "space-time" phase-transition is a first- 
order transition between two dynamical phases which dif- 
fer in their overall "dynamical activity" [HI H] . This tran- 
sition is not controlled by thermodynamic fields. To re- 
veal it one needs to employ a statistical mechanics of tra- 
jectories whereby distributions of trajectories are studied 
by controlling cither dynamical activity or a conjugate 
nonequilibrium field s [81111] . Emergence of dynamic het- 
erogeneity and other fluctuation phenomena can be seen 
as manifestations of this underlying transition [5]. 

In this Letter we show the existence of analogous space- 



time phase transitions in glassy systems driven to non- 
equilibrium stationary states by non-conservative forces. 
In this case the phase-transition is between dynamical 
phases with finite and with vanishing entropy produc- 
tion. We do so by studying ensembles of trajectories via 
the large-deviation method in a class of KCMs driven 
away from equilibrium by non-conservative forces. Such 
forced KCMs [T2J [T3] can be motivated experimentally 
by, for example, pulling a tracer particle through a su- 
percooled liquid or a colloidal suspension P3] . The large- 
deviation method allows to study sub-ensembles of tra- 
jectories with non-typical values of dynamical observ- 
ables in a way that sheds light onto the phase structure 
of the dynamics. We also show that these sub-ensembles 
of rare trajectories can be generated as typical trajecto- 
ries by a mapping to an alternative system with modi- 
fied forces. Such a mapping might allow to bridge the 
conceptual gap between computer generated biased en- 
sembles and the experimental observation of space-time 
phase transitions. 

Biased ensembles of trajectories.- We first recall the 
results for stationary equilibrium dynamics [8j. The ac- 
tivity K [8HTT] is a dynamical (space-time) order pa- 
rameter. Information about the dynamics is encoded in 
the distribution of K, or equivalently, in the dynami- 
cal partition sum Z(s) = (e~ sK )o. The brackets (• ■ -)q 
denote the average over trajectories in the original en- 
semble, i.e., it is a sum over trajectories C(t) weighted 
by P [C(i)]. The field s is conjugate to the activity and 
the factor e~ sK enhances the weight of trajectories with 
either higher than typical activity (s < 0), or with lower 
than typical activity (s > 0). One can think of s as giv- 
ing rise to a biased ensemble of trajectories with weights 
P s [C(t)} oc P [C(i)]e- sX [ c Wl. We refer to this biased en- 
semble as the s-ensemble [5] . 
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Denoting N the system size and £ bs the length of tra- 
jectories, their product Nt b s is the "volume" in space- 
time. For large iVi b s , the dynamical partition sum ac- 
quires a large deviation form |15j with Z(s; Nt b s ) ~ 
e JVtobs^(s)_ The large-deviation function ip(s) generates 
the moments of K in the biased ensemble, e.g., the mean 
activity rate is (K)/(Nt b s ) — —ip'(s) where (• • •) now 
denotes the average in the s-ensemble. The function ip( s ) 
is akin to a free-energy. For KCMs it has a first-order 
singularity at s = indicating a space-time phase tran- 
sition [Sj. This transition is between an ergodic phase 
of finite activity rate and a nonergodic phase of vanish- 
ing activity rate. In the nonergodic phase the system, 
through the kinetic constraints, is able to arrange tra- 
jectories such that the activity grows sublinearly. These 
trajectories "win" at s = + because active trajectories 
with nonzero activity acquire a negative weight in the 
partition sum and thus are exponentially suppressed for 
large i bs- At s = 0~ the opposite occurs and the domi- 
nant phase is the active, ergodic, one. 

We now extend the s-ensemble approach to driven 
systems. Transitions between configurations C occur 
with rates w(C —¥ C) implying the exit rate r(C) = 
^ C iw{C — > C). For a given trajectory C(t) of time 
span i bs the system undergoes K transitions C a _i — > C a 
at times t a where Cq is the initial state, i.e., C(t) = 
(Co, . . . , Crt). The entropy produced in the environment 
in a single transition C — > C is As(C,C) — ln[w(C — > 
C')/w(C — > C)]. In analogy to the activity, we may 
bias trajectories using the time-extensive medium en- 
tropy production s m [C(t)] = ^^ =1 As(C Q _i, C Q ), where 
we sum over all configuration changes. The correspond- 
ing dynamical partition function is 



Z(A;JVt obB ) = < e - A «-)o~e^*W 



(1) 



where A is the parameter conjugate to s m . In analogy 
with the activity, the mean entropy production rate in 
the A-ensemble is given by (s m )/(7Vi b s ) = — W(A). The 
partition sum has a natural transfer matrix represen- 
tation, and \I/(A) is given by the largest eigenvalue of the 
operator [IB] 



W A (C',C) = 



w(C -> C) 
w(C -> C) 



w(C -> C) - r{C)6, 



'CC 



(2) 



In what follows we study the dynamical phase structure 
of driven KCMs by calculating their large-deviation func- 
tions \&(A) from the corresponding operators Wa- 

Driven Fredricks on- Anders en model.- The simplest 
KCM that displays glassy features is the one-spin facili- 
tated Frederickson- Andersen (FA) model [3 [TS] . We in- 
troduce a driven variant of the FA model with evolution 
operator 
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FIG. 1: (a) Mean-field estimate of the large-deviation function 
* (mf) (A) [Eq. M, bold dashed] for the driven FA model with 
c = 0.2, k~/k = 0.7. There is a first-order transition from 
an active (entropy producing) dynamical phase to an inactive 
dynamical phase at A = 0. The "reentrant" transition at 
A = 1 is due to the Gallavotti-Cohen symmetry. The thin 
solid line corresponds to the active branch for all values of A. 
(b) Mean entropy production rate — *'(A). 



For simplicity we consider a bosonic version of the FA 
model, i.e. we allow for multiple occupancy of sites |17j . 
The sum runs over all pairs of nearest neighbors. Here, 
a\ and are creation and annihilation operators at site i 
with number operator n, = ajaj. The first term describes 
the creation and annihilation of excitations at site i with 
rate c = e~^ J , where J is the energy scale and j3 = 
\/{k-QT) is the inverse temperature. The dynamics is 
constrained in the sense that at least one nearest neighbor 
must be excited to allow transitions at site i. The second 
term allows for explicit diffusion of excitations from site i 
to one of its neighboring sites with rates ki^j [19] . Using 
Eq. ([2]), the time-evolution operator for the A-ensemble 
reads 
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We restrict our calculation to a one-dimensional lattice 
with TV sites and periodic boundaries. Further, we con- 
sider spatially homogeneous rates: k + for a hop i — > i + 1 
and kr for i i — 1 such that the steady state solution 
for the mean density becomes independent of the site. 
The largest eigenvalue of ([3| can be estimated in a mean- 
field approximation by maximizing the function W((f>, 4>) 
obtained from the (normally ordered) operator ^ by re- 
placing <v — > $ and a — > <j>. There are two solutions to the 
Euler-Lagrange equations dW/dcf) = and dW/dcf) — 0. 
In terms of the mean density n — <jxj> the two solu- 
tions read: n = and n = (c/16)(3 + y/i + 4k(A)/c) 2 . 
They correspond to the inactive and active phases, re- 
spectively. The A-dependence is contained in 
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FIG. 2: (a) Numerical result for the mean entropy production 
rate in the (2)-TLG for parameters: L — 16, p — 0.5, and 
/ = 0.5. (b) Bimodal distribution of the entropy production 
rate at A = 0.003. The inset shows the kinetic constraint of 
the (2)-TLG: the particle can only move along the arrow if 
the two indicated sites are vacant. 



The mean-field large deviation function becomes 

V ' \ 2n 3 / 2 (Vn -y/c) (A < and A ^ 1) 



(5) 



The function $( mf )(A) and its first derivative are plot- 
ted in Fig. [T] The negative first derivative is the mean 
entropy production rate. It shows a discontinuous transi- 
tions at A = between finite entropy production (A < 0, 
the active phase) and vanishing entropy production (A > 
0, the inactive phase). For A = 1, there is a second "reen- 
trant" transition to an active phase now with a negative 
entropy production, see Fig. [lja). This second transi- 
tion is a reflection of the transition at A = due to the 
Gallavotti-Cohen symmetry, which enforces the relation 
*(1 - A) = #(A) on the large deviation function [T6].l20] 
through Eq. Q . Biasing with A > 1 corresponds to driv- 
ing that essentially changes the sign of all currents due 
to time-reversal symmetry. 

Driven constrained lattice gas. - For a numerical study, 
we turn to a driven variant p~2j US] of a kinetically con- 
strained triangular lattice gas (TLG) [7J [3T] . Lattice gas 
particles diffuse on a two-dimensional triangular lattice 
with periodic boundary conditions. The density of parti- 
cles p is conserved. A particle can only move to a neigh- 
boring site if the two adjacent sites are vacant, see in- 
set of Fig. WCb), This version of the model is called the 
(2)-TLG [2^. In addition, we apply a force / in the x- 
direction. As a consequence the rates for allowed moves 
are w(tp) = e^^ 2 ' Ax ^ where tp = 0, ±7r/3, ±27r/3, 7r is 
the angle between the displacement vector and the x- 
axis. In a single transition a particle moves in the x- 
direction the distance Ax v — 2 cos ip = ±1,±2. We de- 
note by K v [C(t)] the number of moves in the direction 
ip in a trajectory C(t). The activity is then given by 
K[C(t)] = XL -fC^PWL an d the entropy production by 
s m [C(t)} = f± v K v [C(t)]Ax v . 

We use transition path sampling [35] to numerically 
generate biased ensembles of trajectories, for details see 
Ref. [Hj. Fig. [2ja) shows the mean entropy production 



rate (s m ) / (Nt bs) per particle, from ensembles of 10 5 tra- 
jectories for A < 0, and over 10 6 trajectories for A ^ 0. 
There is clear evidence of a transition between a phase of 
finite and a phase of vanishing entropy production sim- 
ilar to Fig. [l|b) for the driven FA model. Compared to 
the mean field solution the transition is smeared out and 
shifted to A > due to finite size effects. The data sug- 
gests that the transition point moves towards A = with 
increasing space-time volume. The probability density of 
the entropy production rate in Fig.[2]jb) for A = 0.003 dis- 
plays a bimodal character, as expected from a first-order 
transition. Just like for the case of equilibrium dynamics 
and the space-time transition in terms of activity [8] , we 
expect the first-order transition between phases of dis- 
tinct entropy production to be generic in driven KCMs. 

Mapping ensembles.- We now address the question of 
whether the biased ensembles of trajectories can be real- 
ized experimentally. In the A-ensemble, the weight of a 
trajectory segment starting in state C, surviving for At, 
and then jumping to state C reads 

We consider long trajectories and include the large devi- 
ation function 'I' (A) for normalization. Collecting terms 
containing At, the exit rate f in a modified dynamics 
must fulfill 

f(C; A) = w{C -> C'\ A) = r(C) + N*(\), (6) 
c 

i.e., the difference between modified and original exit rate 
is given by ^(A) for all states. Simply modifying the 
transition rates as we~ XAs is not enough to give the exit 
rates f. However, adding a sub-extensive term to s m does 
not change the large deviation function. We therefore 
define the modified rates 

W(C -> C; A) = W(C -> C ') e -AAs(C,C')-Au(C,C';A) (?) 

such that the sum of Au along single trajectories is sub- 
extensive in £ b s . E)ue to antisymmetry Au(C',C) — 
—Au(C, C) the sum ^\ Au = vanishes along any closed 
loop £. This property implies O—l independent quanti- 
ties Am, where is the number of configurations C. The 
ft equations ^ then determine {^,Au}. A similar re- 
sult has been found in Ref. [23 for biasing equilibrium 
towards a shear driven ensemble. 

In the simplest case, Au(C,C';X) = u(C';X) — u(C;X) 
with state function u(C; A) depending on A. The sum 
J2 a =i Ait(C a _i,C a ) = u(Ck) — u{Cq) is then a temporal 
boundary term. As we will sec below, u(C) can be inter- 
preted as an energy function. The modified dynamics ^ 
is still Markovian but the rates become nonlocal since a 
transition from site i can now depend on the whole con- 
figuration C. Even for simple KCMs configuration space 
is too big to allow an explicit determination of u(C). 
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FIG. 3: Particle moving in an infinite lattice with alternating 
energy levels (see inset), (a) Symbols are numerical results of 
f (x) / [atoha) for a particle moving with rates Eq. where 
a = — $'(0). The values agree with the normalized mean 
entropy production rate — i S/'(X)/a in the A-ensemble (solid 
lines), (b) The function Aw(A) for / = 2 and c = 0.01. 



To illustrate the idea of mapping the A-ensemble to 
a physical dynamics we consider a simple model related 
to the driven FA model, but without kinetic constraints, 
that can be tackled analytically. A random walker moves 
in a periodic lattice with alternating site energies, see 
inset of Fig. [3|a). The equilibrium rates are wf = 1 
and = c for steps to the right (+) or left (-), where 
c = e~ l3J and J is the energy difference. Applying a 
driving force /, the entropy produced in a single step is 
±/. The total entropy production is s m = fx, where x 
is the distance the particle has traveled. The modified 
rates from Eq. (|7|) are 



w+ = e +//2-A/-A« ) 
,r,- - „-//2+a/-Au 



= ce +//2 ~ A/+A ", 

- „„-//2+A/+A« 



(8) 



with only one independent Au. The two equations ^ 
then lead to quadratic equations for the unknown quan- 
tities Au(A) and \P(A). The solutions are: 



(c - 1) + v/(c-l) 2 +4c[A(A)] 2 



2cA{\) 

*(A) = 2ccosh(//2)L4(A)e Au(A) - 1] 



with A(X) = cosh[/(A - l/2)]/cosh(//2). Again, both 
functions above are symmetric through A(l — A) = A(X). 

Fig. |3][a) shows the mean entropy production in the 
A-ensemble of the original driven system with rates Wi 2 - 
Around A — and A = 1 we observe a crossover be- 
tween regimes of high entropy production and one of 
low entropy production (something similar has been ob- 
served for a driven particle moving in a periodic poten- 
tial [21]). Fig. |3ja) also shows that the biased entropy 
production coincides with the typical entropy production 
in the mapped problem with rates wf 2 from Eq. (j^J). 
The driving force in the mapped system transforms as 
/ = /(l — 2A). The barrier height is adjusted through 
Au [see Fig. [3|b)] and the barrier becomes maximal for 
A = 1/2 and vanishes for large |A| with Am — > f3J/2 to 
allow for the maximal current. One should note that the 



correct expression to compare — ^'(A) with is f(x) [where 
the brackets are now the average over the dynamics with 
rates ^] and not the actual entropy production f(x) of 
the mapped system. 

Conclusions.- We have shown that the first-order 
space-time phase transitions between active and inactive 
phases observed in the equilibrium dynamics of model 
glasses [8] are also present in systems driven to non- 
equilibrium stationary states. In this case the transition 
is between phases with differing entropy production rates. 
A next step is to search for evidence of such phase tran- 
sitions in driven atomistic liquids in analogy with what 
was done in Ref. [5] for the case of equilibrium fluids. 
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